Spin chirality induced by the Dzyaloshinskii-Moriya interaction 
and the polarized neutron scattering 
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We discuss the influence of the Dzyaloshinskii-Moriya (DM) interaction in the Heizenberg spin 
chain model for the observables in the polarized neutron scattering experiments. We show that differ- 
ent choices of the parameters of DM interaction may leave the spectrum of the problem unchanged, 
while the observable spin-spin correlation functions may differ qualitatively. Particularly, for the 
uniform DM interaction one has the incommensurate fluctuations and polarization-dependent neu- 
tron scattering in the paramagnetic phase. We sketch the possible generalization of our treatment 
to higher dimensions. 



Since the works by Dzyaloshinskii and MoriyaB, the 
antisymmetric spin exchange interaction plays an im- 
portant role in the physics of condensed matter. Be- 
ing introduced for the explanation of the weak ferro- 
magnetism in antiferromagnets without center of in- 
version, the Dzyaloshinskii-Moriya (DM) interaction is 
found nowadays in various problems of magnetism and 
statistical physics. 

Being the relativistic effect, the magnitude of DM in- 
teraction, D, is generally expected to be small in coat 
parison with the usual symmetric superexchange, J Era 
In some compounds, however, this interaction can attain 
a sizeable value. For instance, one. has D/J = 0.18 in 
the hexagonal pexovskite CsCuCi3tJ and D/Jm 0.05 in 
copper benzoate.cl 

Remarkably, the DM interaction in the latter com- 
pounds takes place in the quasi-one-dimensional spin sub- 
systems. On this reason, we are primarily concerned be- 
low with the one-dimensional (ID) situation of a quan- 
tum spin chain. 

Generally, the DM interaction between two spins, Si,2j 
is written as D(Si x S2) with an axial DM vector D. In a 
chain, D may spatially vary both in direction and mag- 
nitude, however, the symmetry arguments usually rule 
out most of the possibilities and confine the theoretical 
discussion to two principal cases. The first one is the 
uniform DM interaction, D = const over the syatem.EI 
The second case is the staggered DM interaction,!!! with 
antiparallel D on adjacent bonds. 

Among the other studies, we should mention the dis- 
cussion of the XY spin chains with randomly distributed 
values of Da and the growth models with imaginary uni- 
form D = iX leading to non-Hermitian Hamiltonian.tj A 
model of XY spin chain with a ternary DM interaction 
was introduced and solved recently.Q 

In a present paper, we deal mostly with two above 
cases, uniform and staggered DM interaction. We con- 
sider also a model of a non-ideal lattice, where one finds, 
say, an almost uniform situation with one possible DM 
value, D, taking place on a chain fragment of average 
length li, and another value, — D, on a chain fragment of 
length I2, while Z1.2 3> 1. The situation is then described 
in probabilistic terms. We show that, being the situa- 



tion uniform, staggered or random, the spectrum of the 
Heisenberg chain with the DM interaction is equivalent 
to one of XXZ spin model and is computed exactly for 
spin 1/2. Of course the observable susceptibilities can 
differ crucially, as we demonstrate below. 

Therefore we extend the previous result by Alcaraz and 
Wreszinski, that the Heisenberg ID model with the uni- 
form DM interaction is reduced to XXZ spin exchange 
model and is exactly solvableE 

For the uniform and almost uniform DM interaction, 
we show that the observable spin-spin correlation func- 
tion possess an incommensurate structure. This in- 
commensurability phenomenon was noted previously for 
the XY spin chain model and uniform antisymmetric 
spin interactions We show further that in this case the 
spin susceptibility tensor x Q/3 acquires an antisymmetric 
part. This leads to the appearance of the polarization- 
dependent part of the neutron scattering cross-section, 
which makes possible the direct observation of the direc- 
tion and the value of the DM vector D. Our results are 
applicable in the absence of the long-range magnetic or- 
der in the system. They also can be generalized towards 
higher-dimensional situation, as discussed below. There- 
fore our treatment may provide an explanation of the ear- 
lier experiments in cubic ferromagnet MnSi,E3 where the 
DM-induced incommensurability of the magnetic fluctu- 
ations and polarization dependence of the neutron scat- 
tering were observed both below and above the Curie 
ordering temperature. 

We consider the spin chain Hamiltonian of the form 



H = J2(JS l S l+1 +T> l [S l x S m ]) 



(1) 



1=1 



with AF Heisenberg coupling J and the Dzyaloshinskii- 
Moriya term D/ . We choose the vector D/ to be directed 
along the z— axis. 

We observe that Ti is simplified upon a canonical trans- 
formation TL 



He lU with 



U = ^aiSf, at =J2 tan_1 (A/ J), 



(2) 



1=1 



1 



and ai — 0. Note that this transformation works for 
all values of 5. The periodic boundary conditions (BC) 
require ul+\ = 0mod27r which relation is not gener- 
ally satisfied. However, in the thermodynamic limit 
L — > oo the influence of BC can be negliged.B Introducing 
Sf = S x ± iS V j, one can easily see that 



S± = e- iU S^e iU = S±e* ia < , Sf = S 



i ■ 



(3) 



and our choice of the coefficients ai removes the antisym- 
metric part of the Hamiltonian, [S; x S;+i]. We consider 
below two principal possibilities: the "uniform" situation 
D[ = JtanJ and the "staggered" one D\ — (— 1)' JtanJ. 
In both cases the Hamiltonian is reduced to the XXZ 
model : 



H = y{J x {SfSf +l + SfSf +1 ) + J5f 5f +1 ) 



1=1 



(4) 



with J x = y/ J 2 + Df independent of I. 

Our subsequent discussion is based on the observa- 
tion, that the DM interaction results in two effects for 
the observable susceptibility of the system. First ef- 
fect is the modification of the spectrum, as seen in 
the equivalent Hamiltonian The appearance of the 
"easy-plane" anisotropy (J x > J), however, does not 
lead to a gap in the spectrum. The exact solution of 
(Q) for 5 = 1/2 showstH that the correlation functions 
(5f5* ) ~ |Z - m|-" and (S?S* m ) ~ |Z - m^ 1 '" with 
v = 1 — \S\/tt. Since the value of DM exchange is ex- 
pected to be small, \Di\ -C J, the long-distance decay 
of the above correlation functions is described to a good 
accuracy by the "isotropic" Heisenberg situation, with 
v = 1. The second effect of DM interaction for the ob- 
servables is the explicit dependence of the relation (|^) 
between the new and old spin variables on the values of 
Di. We focus our attention below on the latter effect, 
which leads to the qualitative changes in the experimen- 
tally observable susceptibilities. 

The two-time Green's function for the operators A and 
B is defined as XAEi(t) = -i0(t){[A(t), B\) where [...,...] 
stands for a commutator and 9{t) = 1 at t > 0. 

Upon the "twist" e lU the ^-component of spin oper- 
ators remains unchanged, and one has for the longitu- 
dinal zz susceptibility xj£(t) = -iO(t){[Sf (*), S 7 m ]) = 

-i6{t){[Sf{t),S^\) = gf m (t). Therefore the observable 
Xfm(t) nas a commensurate antiferromagnetic modula- 
tion. 

The expressions for the transverse spin susceptibility 
are mqse complicated. It is convenient to introduce the 
matrixEa 



xL(t) = -i9(t) 



([5f(*),5*]>, ([5f(i),5*]) 
([5?(t),5*]>, ([Sf(t),sy m ]) 



(5) 



in the initial system (Q). In the simpler "twisted" system 
(i we have -iQ{t){[Sf (*), &]) = -W)([Sf (*), &J) = 



GUt) and ([Sf(t),S*]) = ([5f(i),5*]) = 0. Returning 
back to quantities {[Sf {i), 5^J) with the use of (0), we 
get 



xL(t) = QL(t) 



cosaz,™, -sinaj. 
smai, rn , cosai„ 



(6) 



with a ijm = a; - a m . 

For later comparison, it is worth to consider 
first the case of the staggered DM interaction.Q We 
have Di = {—l) l Jt&n8 and ai, m = ((—1)' — 
(— l) m )<5/2. In this case we write cosa^m = 
cos 2 (<5/2) + (-l) ( - m sin 2 (5/2) and sina; >m = ((-1)' - 
(— l) m )(sin(5)/2. Clearly, the off-diagonal components 
xfm(*) °^ ^he ma t r i x © d° n °t depend on the differ- 
ence (Z — m) only and the two-momenta Fourier trans- 
form A(q,q',u) = J dtJ2i m e lql ~ lq ' m ~ liAJt Ai m (t) should 
be introduced. Then we obtain the off-diagonal com- 
ponents in the form x cy {q, Q'i w ) — (G J ~{q,^j) ~ G {q — 
it, lj)) J2 T 6{q — q' — n + t) where J2 T stands for the sum 
over all vectors r = 27m of the reciprocal lattice. How- 
ever, apparently in all physical observables one finds the 
symmetrized form of the susceptibility {q — q') and the 
off-diagonal terms in the matrix x ± vanish. Therefore in 
the case of staggered DM interaction one is left with the 
diagonal component of the matrix X fe^) °f the form : 



X X {q,u) = G ± {q,Lo) cos 2 {5/2) 

+t;[S ± (q - tt, w) + Q X {q + n, u)] sin 2 (<5/2) 



(7) 



We see that the regions of the AF and the ferromag- 
netic fluctuations are mixed in the observable suscepti- 
bility. A consequence of this featureo is the anomalous 
temperature behavior of the uniform static susceptibility 
X xx {0,0) = x vv {0,0) for the AF chain (see also the dis- 
cussion aftex_Eq.(5.3) in the original Moriya's papera). 
It is knowrcj that in the Heisenberg 5=1/2 chain one 
has £^(0,0) - J- 1 and G ± {ir,0) ~ T- 2+l/ , therefore 
X zz (0,0) - J- 1 and 



: (0,0)= X ra (0,0) 



J- 1 [const + S\J/T) 



1+|<5|/7T1 



(8) 



The Eq. (g) has a simple physical meaning. Indeed, in 
the considered case the operator U "cants" the local co- 
ordinate frames by an angle ±5/2. It leads effectively to 
the non-compensated spin AS = 5/4 in the x — y plane. 
Being the spins AS free, it would then lead to the Curie 
law for the susceptibility x ~ {AS 2 )/T. The ID charac- 
ter of the interacting spin system results in the nontrunal 
exponent in the T— dependence of this term (cf. alsoa). 

At the same time, the above transfer of the spectral 
weight, Eq. [j], is apparently negligible to be observed in 
the neutron scattering experiments. 

On the other hand, for the "uniform" DM interaction 
Di = D we have aj, m = (Z — m)5. It results in the incom- 
mensurability of the transverse spin correlations. Fourier 
transforming Eq. (o), we obtain 



2 



5,cj) 



1, 



1 



(9) 



Let us discuss the physical consequences of this expres- 
sion. Evidently, ~x xx (Q,Q) ~ const in this case and the 
presence of the "uniform" DM interaction is not revealed 
by the measurements of the temperature dependence of 
the uniform static susceptibility. Much more interest- 
ing are the implications of (^|) for the polarized neutron 
scattering experiments. The basic quantity here is neu- 
tron scattering cross-section, which is connected to the 
Green's function x a ^(qT UJ ) of the spin system. It is con- 



venient to write x a ^ = Xs P ~ *Xa > w hh the symmetric 
and antisymmetric tensors, Xs^ an< ^ fcA ^ respectively. 
Up to fundamental constants, we haveEf : 



(10) 



where N(ui) is the Planck function, the unit vector 
Q = Q/Q is directed along the neutron's momentum 
transfer Q, e Q( a 7 is totally antisymmetric tensor, Po is 
the incident neutron's polarization and q is the on-chain 
projection of Q. 

From (10) we see that if the whole crystal is character- 
ized by Dzyaloshinskii vector D (uniform situation, Eq. 

then the polarization-dependent part of cross-section 
is non-zero and is given by 



dQdiu 



(DQ) (QP ) Im 



g ± (g + S,cj)-g ± (q-5,cu) 
2D 



(11) 



A certain subtlety should be discussed here. Under 
the parity transformation we have q —* — q, D — > D, 
Po — > Po- At the first glance S — ► S and thus d 2 a\ / dVldio 
changes the sign, as it should not be. An inspection of 
@ , (0) (^|) shows however that the quantity 8 in (^lj) ap- 
pears as the differential of ai . The latter object is the sum 
of the phases Sj over the bonds j to the left of /. Hence we 
have under the parity transformation cti — ► — ai + const 
and 5 — > —S in (p"l|), which restores the desired property 
of the cross-section. 



The contribution of the symmetric part of x 1° 
d 2 a(q,uj)/dflduj in the considered case of uniform D is 
two-fold. One still has the commensurate fluctuations 
of spin components S^ljD, with a peak at the AF posi- 
tion. At the same time, the DM interaction splits the AF 
peak related to transverse fluctuations into two peaks of 
the weight 1/2, Eq. |^. The relative weights of these two 
structures depend on the direction of Q as follows 

<*Wg,") „ (1 + q 2) Im GHq + 6,u) + GH<i-8,w) 

dfldu> " 2 

+{l-Ql)ImQ\\q,u). (12) 



Note that at low temperatures, the incommensurate 
peaks have more singular behavior according to our dis- 
cussion after Eq. (||). 

An important thing to be stressed here is the follow- 
ing. It is known that the incommensurate long-range 
magnetic structures may arise due to the competing in- 
teractions in the spin system. In this case one expects 
that all three diagonal components of the spin suscep- 
tibility x aa are peaked in the paramagnetic region at 
the same incommensurate wave-vector. The off-diagonal 
components of x a ^ are absent. This is fairly different 
from the picture described above, Eqs. [ll], [l^. Hence the 
experimental observation of the incommensurability phe- 
nomenon in the paramagnetic phase, accompanied by the 
polarization dependence of the neutron scattering cross- 
section could serve as an indication to the presence of 
the uniform DM interaction. Remarkably, the value and 
the direction of the pseudo- vector D can be, in principle, 
determined this way. 

In reality, however, the macroscopic sample is rarely 
uniform and it should be expected to split to domains 
with different directions of D. To account for this sit- 
uation, it is instructive to analyze a model where the 
value of the DM interaction Di takes randomly two val- 
ues ±Jtan(5. 

Consider first the oversimplified case when ((-D/)) = 
and {(DiD m )) = for I ^ m, here ((. . .)) denotes aver- 
aging over the realizations. In this case the spectrum 
is still defined by Eq.(^), and x zz is given by the above 
expression. At the same time, one can easily show that 
the averaged susceptibility x 1 ' has a diagonal form and 
exhibits an exponential decay of correlations : 



((xL(t))) = GL(t)eM-\l-m\/h 



(13) 



with the correlation length = — l/ln(cos5) ~ S~ 2 . 

Now consider a more realistic situation when one still 
has Sj = ±5, but the signs of Sj on the adjacent bonds 
are correlated. The diagonal and off-diagonal parts of 
the matrix (^J) are given, respectively, by the real and 
imaginary part of the average 



i-i 



,(exp ia it7 



,S L )exp(i ^2 h)- 



(14) 



k—m 



We assume that the joint distribution function 
p{5\, . . . , S n ) has a Markovian character, p(Si, . . . , 5 n ) — 
p(S\, . . . , S n -i)p(5 n \S n -i). In this physically important 
case we arrive at the diekeitomous Markovian noise <5i 



with a discrete "time" We set ((Sj)) — 5d which 
defines the on-site ("equilibrium") probability as po — 
(^2> ^T/" The matrix p(S n \S n -i) satisfies the "conser- 
vation laws" for the total and equilibrium probabilities, 
(1,1) ■ p(S n \S n -i) = (1,1) and p(S n \S n -i) ■ p = Po, re- 
spectively. These equalities fix p(S n \5 n -i) in the form 



l-ar(l-d), x(l + d) 
x(l-d), l-x(l + d) 



3 



for all n, which corresponds to the following correla- 



tor on the adjacent sites 



S 2 (l — d 2 )(l — 2x). The latter equalities mean that the 
absence of correlations corresponds to x — 1/2 and the 
correlation lengths for positive and negative sequences 
of Sj are I/Z12 — ^(1 T d). Introducing the matrix 
T> = diag(e lS , e~ lS ), the quantity (|l4|) is represented as a 
product (1, 1) ■ (T> ■ p)'-™ 1 - 1 -T>-p , which is evaluated us- 
ing the multiplication rules for the Pauli matrices. After 
straightforward, though tedious, calculation, we obtain 
the average ( |l4| ) in general form, which is somewhat sim- 
plified in two principal cases of small 5 : i) 5 <C x ~ 1 
and ii) 5 ~ x <C 1. In the first case, keeping the terms of 
order of <5 2 , we have 

((e m "*+"' m )) ~ exp[n(iSd - 6 2 ai ) + S 2 a 2 ] (15) 
- <5 2 a 2 (l - 2x)"exp[(n- l)(~iSd + S 2 ai )} 

here ai = (1 - x)(l - d 2 )/(2x) and a 2 = (1 - d 2 )(l - 
2x) I (4a: 2 ) . We see that the incommensurability wave vec- 
tor is defined by the average on-bond value 8d. Note that 
Eqs. (m), ( |l3| ) are recovered at d = 0, x = 1 and d = 0, 
a; = 1/2, respectively. 

When x ~ S -C 1, we come to a more complicated 
situation. We have in this case 



, , sinh bn . . . n 
cosh on -I — (x + tod) 



(16) 



with b = y/x 2 + 2ix8d- S 2 . We return to Eq. @ at 
x = and d = 1. 

It should be stressed that at d = one has 
/ m (( e <om+»,m)) = both in Eqs. ©, (jl|) and in gen- 
eral case. It corresponds to the fact that the off-diagonal 
components of susceptibility x^> Eq. (JsJ) , vanish in the 
system with zero average Dzyaloshinskii vector ((D)). 
As a result, the observable cross-section is polarization- 
independent. The position of maximum of the transverse 
spin fluctuations, though, may be incommensurate one 
for the almost uniform DM interaction, as seen from Jl6| ) 
at d = and x — > 0. 

Now we discuss the possible generalization of our ap- 
proach to a higher dimensional case. Consider a planar 
system with spins S; m labeled by two indices. The inter- 
action between spins takes place in two directions, and we 
write the corresponding quantities as and , with 
a = x, y. For simplicity we consider the case when the 

vectors lie along one direction, with possible varia- 
tion in their magnitude. We introduce then two angles, 



I in 



tan-l(D«/jfi?) and 5, 



W - 



tan-H^/J^). 



We are interested to arrive to a symmetrized Hamil- 
tonian, similar to (Q), by making the transformation 
with U = ^2^ m =i&imS[ m . One can show that this 
transformation is possible if and only if Sf. 



l,m+l ~~ °lm 



Sf +1 m — Sf m . In this case oti m is uniquely defined and 

may be written in a form ai m = $^=1 +Y<T=i ' 
The latter relations are not surprizing when we note that 



in the continuum limit they read as V x S r = 0, while 
5 r = ((5^,(5^) = Va r . In particular, the two-dimensional 
situation with D 



(x) 
Irn 



(v) 

D and D\^[ — allows this trans- 



formation. The calculation of observables and the gener- 
alization to a three-dimensional case are done in the way 
similar to the above one. 

In conclusion, we discuss the observables in the spin 
system with the antisymmetric DM interaction. We show 
that in one spatial dimension the exactly found spectrum 
of such problem may coincide for different choices of the 
parameters of DM interaction. Despite this fact, the ob- 
servable spin-spin correlation functions may differ cru- 
cially, and we discuss this feature with the application 
to the neutron scattering experiments. In particular, the 
incommensurability and the polarization dependence of 
the neutron scattering may be used for the determination 
of the value of the uniform DM interaction. 
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